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ABSTRACT 


An  investigation  of  power-density  autospectrnni  estimation  by  means  of  overlapped  Fast  Fourier 
Transform  (FFT)  processing  of  windowed  data  is  conducted  for  four  candidate  spectral  windows  with 
good  side-lobe  behavior.  A  comparison  of  the  four  spectral  windows  is  made  on  the  basis  of  equal 
half-power  resolution  bandwidths.  The  criteria  for  comparison  are:  (1)  statistical  stability  of  the 
spectral  estimates,  (2)  leakage  (side  lobes)  of  the  spectral  windows,  (3)  number  of  FFTs  (number  of 
overlapped  pieces)  required,  and  (4)  sizeof  each  FFT  required.  The  dependence  of  these  criteria  on 
the  amount  of  overlap  is  investigated  quantitatively. 

Some  striking  invariances  are  discovered.  Specifically,  it  is  shown  that  the  ultimate  variance- 
reduction  capabilities  of  the  four  windows,  as  measured  by  the  equivalent  number  of  degrees  of 
freedom  (EDF),  are  virtually  identical  under  the  constraint  of  equal  half-power  bandwidths.  Further¬ 
more,  when  the  proper  overlap  is  used  for  each  window,  the  stability  of  this  method  of  spectral  esti¬ 
mation  is  identical  to  that  of  the  "indirect"  correlation  approach.  Also,  the  number  of  FFTs  required 
to  realize  99  percent  (or  less)  of  the  maximum  EDF  is  virtually  independent  of  the  particular  window 
employed.  The  required  fractional  overlap  of  the  four  data  windows  lor  99  percent  (or  less)  of  the 
maximum  EDF  is  virtually  independent  of  the  product  of  the  available  time  and  the  resolution  band¬ 
width,  although  it  dots  depend  on  the  particular  window.  Tables  of  required  overlap  are  presented. 
The  only  tradeoff  among  the  four  windows  is  that  those  withbetter  side  lobes  require  larger-size  FFTs. 
All  of  these  results  are  derived  for  a  Gaussian  random  process,  under  the  assumption  that  the  resol¬ 
ution  bandwidth  of  the  spectral  window  is  smaller  than  the  finest  detail  in  the  true  spectrum. 

Pules  of  thumb  for  the  maximum  EDF  and  the  number  of  FFTs  required  tv'  realize  99  percent  of 
the  maximum  EDI  are  given.  The  possibility  of  weighting  individual  spectral  estimates  unequally 
in  order  to  optimize  the  EDF  is  investigated;  the  gain  is  found  to  be  negligible  for  cases  of  practical 
interest. 
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An  investigation  of  power-density  autospectrum  estimation  by  means  of  overlapped  Fast 
Fourier  Transform  (FFT)  processing  of  windowed  data  is  conducted  for  four  candidate  spectral 
windows  with  good  side-lobe  behavior.  A  comparison  of  the  four  spectral  windows  is  made  on 
the  basis  of  equal  half-power  resolution  bandwidths.  The  criteria  for  comparison  are:  (1)  sta¬ 
tistical  stability  of  the  spectral  estimates,  (2)  leakage  (side  lobes)  of  the  spectral  windows, 

(3 1  number  of  FFTs  (number  of  overlapped  pieces)  required,  and  (4)  size  of  each  FFT  required. 

Some  striking  invariances  are  discovered  in  the  ultimate  variance-reduction  capabilities 
and  in  the  number  of  FFTs  required  to  realize  99  percent  (or  less)  of  the  maximum  EDF.  The 
onL  traclc- off  among  the  four  windows  is  that  those  with  better  side  lobes  require  larger-size 
FFTs. 
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SPECTRA I.  ESTIMATION  HY  MEANS  OF  OVERLAPPED 
FAST  FOURIER  TRANSFORM  PROCESSING  OF  WINDOWED  DATA 


INTRODUCTION 

Estimation  of  the  power-density  spectra  of  stationary  random  processes  is 
an  important  problem  and  occurs  frequently  in  many  fields.  The  resolution  of 
closely  spaced  frequency  components,  with  limited  amounts  of  data,  presents 
inherent  limitations  on  the  statistical  stability  of  the  estimates.  Also,  the  pre¬ 
vention  of  leakage  of  undesired  frequency  components  into  the  frequency  being 
analyzed  dictates  a  careful  choice  of  data  weighting.  Lastly,  the  extent  and 
complexity  of  the  data  processing  required  to  realize  the  desired  resolution, 
stability,  and  leakage  control  are  important  considerations. 

The  fundamental,  conflicting  desires  involved  in  spectral  estimation  be¬ 
come  painfully  obvious  when  the  amount  of  data  available  for  analysis  is  limited 
and  can  not  be  augmented  by  additional  measurements.  For  example,  the  avail¬ 
able  record  length  may  be  limited  by 

a.  nonstation  ary  conditions  (changing  environment), 

b.  storage  limitations, 

c.  equipment  failure,  and 

d.  time-sharing  requirements. 

Although  factors  b.  c,  and  d  can  often  be  remedied  or  corrected,  factor  .. 
often  can  not  be  controlled.  Thus,  only  a  small  segment  of  the  time  record 
may  be  usable  for  caeh  spectral  analysis.  If  fine  frequency  resolution  is  de¬ 
sired,  the  limited  number  of  independent  observations  available  makes  stable- 
estimation  impossible  in  sonic  cases.  One  must  then  be  willing  to  accept 
coarser,  but  more  stable,  spectral  estimates. 

The  two  fundamental  parameters  that  control  the  performance  of  spectral 
estimation  are  the  available  record  length,  T,  in  which  the  sample  of  the 
random  process  is  assumed  stationary,  and  the  desired  frequency  resolution, 

B,  of  the  spectral  analysis.  Large  values  of  the  fundamental  BT  product 
yield  good  performance  of  the  analysis  teclinique,  but  small  values  arc  often 
forced  upon  us  by  too  small  a  record  length  T  or  too  fine  a  desired  resolution 
3.  The  problem  here  is  to  make  maximum  use  of  the  available  data. 
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Spectral  analysis  has  received  much  attention  in  the  past  [1-6],  especially 
since  the  advent  of  the  Fast  Fourier  Transform  (FFT)  [7,8].  In  particular,  the 
method  of  averaging  short  modified  periodograms  [9]  is  a  prime  candidate  for 
spectral  analysis  —  for  several  reasons.  First  of  all,  nonstationary  trends  in 
the  data  are  more  readily  observable  through  the  time-local  spectral  estimates 
of  each  segment.  Second,  the  size  of  each  FFT  can  be  kept  reasonably  small, 
thereby  reducing  storage,  execution  time,  and  round-off  error.  Third,  the  fre¬ 
quency  resolution  is  easily  controlled  by  the  choice  of  segment  length,  and  leak¬ 
age  (side  lobes)  can  be  controlled  by  the  proper  choice  of  window  in  the  time 
domain.  Lastly,  overlapped  segments  of  windowed  data  utilize  more  fully  the 
variance-reduction  capability  of  a  given  record  length. 

The  problem  to  be  addressed  here  has  to  do  with  the  choice  of  window  and 
amount  of  overlap  to  employ  for  a  particular  application.  Specifically,  if  we 
employ  a  window  with  very  small  spectral  side  lobes,  how  much  should  the 
segments  be  overlapped,  and  does  the  overlap  vary  greatly  with  the  particular 
window  selected?  (Fifty-percent  overlap  has  been  suggested  as  a  reasonable 
procedure  for  the  triangular  data  window  [9,  p.  72].)  How  many  FFTs  of  what 
size  have  to  be  performed  for  the  different  windows?  Is  the  variance-reduction 
capability  dependent  on  the  particular  window? 

Four  windows  will  be  investigated.  They  are  called  data  windows  in  the 
time  domain,  where  they  are  multiplied  by  the  available  data  record;  they  will 
be  called  spectral  windows  in  the  frequency  domain,  where  their  main  effect 
enters  via  convolution.  The  four  data  windows  are  called  triangular,  cosine,  * 
quadratic,  and  cubic  and  will  be  documented  in  a  later  section.  Hamming 
weighting  [2,  p.  14],  although  it  possesses  good  adjacent  side  lobes,  is  not 
considered  here  because  the  spectral  window  decays  very  slowly  with  frequency, 
thereby  responding  to  frequencies  far  removed  from  the  analysis  band  of 
interest. 


PROBLEM  DEFINITION 

Consider  that  stationary  random  process  x  has  been  observed  for  a  time 
interval  T  seconds;  that  is,  x(t)  for  0<t<T  is  available.  Let  the  power- 
density  autospectrum  of  this  process  at  frequency  f  bommnted  by  G(f),  where 
double-sided  spectral  notation  will  be  employed.  Wc  wish  to  estimate  spec¬ 
trum  G  with  a  resolution  of  13  hertz,  where  B  is  the  half-power  (-3  dB) 
bandwidth  of  the  desired  resolution. 


'Also  called  Hanning  weighting. 


DAT  A  WINDOW  w 


The  method  of  obtaining  the  spectral  estimates  is  depicted  in  Fig.  1.  A 
data  window  w  of  duration  L  seconds  is  applied  successively  to  the  available 
data  x  in  the  overlapping  intervals  (0,L),  (S.S+L),  ...  ,  ((P-l)S,  (P-1)S+L). 

S  is  the  amount  of  shift  each  adjacent  data  window  undergoes,  and  P  is  the 
total  number  of  pieces  or  segments  employed.  Since  only  T  seconds  of  data 
are  available,  we  must  have 


(P-1)S+L  <  T  .  (1) 

The  segment  length  L  should  be  large  enough  that  the  correlation  function  of 
process  x  is  effectively  zero  for  delays  larger  than  L/2.  (The  relation  be¬ 
tween  frequency  resolution  B  and  segment  length  L  is  discussed  quantitatively 
later.)  The  form  of  the  data  window,  depicted  in  Fig.  2,  is  even  about  the 
origin  and  real.  Also,  w(t)  is  zero  for  iti>L/2. 
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Fig.  2.  Data  Window 


When  the  overlap  in  Fig.  1  is  a  significant  fraction  of  the  segment  length, 
the  effective  use  of  the  available  data  x  is  fairly  uniform  over  the  entire  inter¬ 
val  T  except  for  the  edges  of  the  data,  where  a  gradual  taper  over  an  interval 
of  length  L/2  takes  place.  This  is  consistent  with  earlier  suggestions  [lO, 
p.  58]  for  maximum  use  of  the  available  data.  Notice  that  the  percentage  of 
taper  depends  on  the  desired  frequency  resolution  and  available  record  length, 
and  is  not  a  constant,  such  as  10  percent,  as  lias  occasionally  been  suggested. 

The  estimate  of  the  power-density  spectrum  G  is  obtained  as  follows. 
First,  a  Fourier  transform  on  the  p-th  windowed  section  is  performed*: 

Y^(f)  =  J“dl  exp(-i2jrft)  x(t)  wj\  -  j-  -  (p- l)sj  ,  1  <  p  g  P  .  (2) 

The  spectral  estimate  at  frequency  f  is  then  available  as  the  average  of  the  P 
pieces: 


Equation  (2)  assumes  a  continuous,  rather  than  discrete,  form  of  signal 
processing  [2,  Secs.  4-11  versus  12-2 1] .  However,  if  the  discrete  version  of 
Eq.  (2),  where  samples  of  x  are  taken  At  seconds  apart,  is  such  that  aliasing 
is  negligible,  there  is  little  difference  between  the  two  methods  of  spectral 
analysis  [2,  pp.  37-39  and  123-125].  We  shall  assume  that  At  is  so  chosen 
and  confine  attention  here  iO  the  continuous  processing  technique  of  Eq.  (2).  Of 
course,  in  practice,  Eq.  (2)  is  approximated  by  a  discrete  Fourier  transform 
[9,  p.  70]  ,  in  which  case  dc  and  linear-trend  removal  should  he  considered  for 
the  sampled  data  [2,  pp.  47-49], 

The  spectral  estimate  ^(f)  in  Eq.  (3)  is  a  random  variable.  Its  mean  and 
variance  are  evaluated  in  Appendix  A  under  the  assumptions  that  x  is  a  Gauss¬ 
ian  random  process  and  that  the  frequency  resolution  of  the  spectral  window 
|  W|  “ ,  whe  re 

W(f)  =  /dt  cxp(— i2 a  ft )  w(t)  ,  (4) 

is  narrower  than  the  finest  detail  in  the  true  spectrum  G.  (This  latter  assump¬ 
tion  is  equivalent  to  that  given  under  Eq.  <  1  >  for  segment  length  L.)  The  re¬ 
sults  are 

•Integrals  without  li:  1  its  are  over  the  range  of  non-zero  integrand. 
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(3) 


El 6(f)}  ,/d.  G(f- .')  |W(»)|2SG(f) / d*  |W(^)12  , 

Var  f  G(f)}  S  ^  £  ^1  -  /m  exp(i2rukS)  G(m)  |\V(f-*«)|S 


P-1 

E 

k=-(P-l) 

P-1 

E 

k=-(P-l) 


(6) 


where 


<t> ^  (r)  = dt  w(t>  w*(t  -  r)  . 


(7) 


Relation  (5)  shows  that  the  mean  of  the  spectral  estimate  is  equal  to  the  con¬ 
volution  of  the  true  spectrum  G  with  ihe  spectral  window  |W|2.  Relation  (6) 
expresses  the  variance  of  the  spectral  estimate  in  terms  of  the  number  of 
pieces  P,  the  shift  S,  and  the  autocorrelation  <t>w  of  the  data  window .  The 
result,  Eq.  (G),  holds  if  f  is  greater  than  the  bandwidth  of  the  spectral  win¬ 
dow,  the  right  side  of  Eq.  (G)  must  be  doubled  if  f  =  0  [see  also  9,  p.  71]  . 


The  equivalent  number  ot  degrees  of  freedom  (EDF)  in  spectral  estimate 


G  is  defined  as  [2,  p.  22] 


E2  {6(f)} 
Var  {6(f)} 


K  =  2 


2P 


P-1  /  x 

E  (‘  -  t) 


k--(P-l) 


$H.(kS) 


V°) 


2  ’ 


(3) 


employing  Eqs.  (4)  through  (G).  Notice  that  under  the  assumptions  given  above, 
K  is  independent  of  the  value  of  frequency  f  and  true  spectrum  G;  for  f  =  0, 
K  is  given  by  one  half  of  Eq.  (8). 

For  computational  purposes,  it  is  convenient  to  define  a  normalized  data 
window  u  according  to 


(9) 


:> 


u(x)  -  w(Lx)  . 


Then  u  contains  the  shape  Information  of  data  window  w  but  extends  only  over 
the  interval  (-1/2,  1/2).  It  follows  that 


0  (r)  =  L  <f>  (t/L)  , 
w  u 


where 


Vr,“/‘ 


dt  U(t)  U*(t-  r)  . 


The  EDF  in  Eq.  (8)  then  becomes 


P-1  /  4>  (k 

Z  1  T 

t  — fP-1)  U 


<t>u(kS/L) 
~  (Of 


In  order  to  minimize  the  fluctuations  in  the  spectral  estimate  G,  we 
should  maximize  I<  in  Eq.  (12).  To  accomplish  this  for  a  given  number  of 
pieces  P  and  segment  length  L  the  shift  S  in  Eq.  (12)  should  be  chosen  as 
large  as  possible  so  that  4>u  is  a;.'  small  as  possible  at  S/L,  2S/L,  etc. 
However,  since  Eq.  (1)  dictates  a  constraint  among  these  variables,  the  best 
choice  of  shift  S  —  for  given  record  length  T,  number  of  pieces  P,  and 
segment  length  L  —  is  given  by  equality  in  Eq.  (1): 


S  -  (for  P  >  2)  . 


Substituting  Eq.  (13)  in  Eq.  (12),  we  have 


/ i  a 

7 

P-1  /  ,  \ 

*  kP  J 

(i-~) 

u  \  P  -  1/ 

L*  V  p/ 

<p  (0) 

k=-(P-l) 

u 

Equation  (14)  expresses  the  EDF  as  a  function  of 
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P,  number  of  pieces  in  the  average, 

T/L,  ratio  of  record  length  to  segment  length,  and 

<t>  ,  autocorrelation  of  shape  ol  data  window 
u 

The  problem  now  is  to  maximize  the  EDF  in  Eq.  (14)  by  choosing,  subject  to 
specified  recced  length  T  and  desired  half-power  frequency  resolution  B,  the 
number  of  pieces,  P,  for  several  data  windows,  w.  It  will  turn  out  that  the 
optimum  value  for  P  is  not  infinite. 

One  special  case  of  Eq.  (14)  is  worth  comment:  if  P  <  T/L,  the  values 
of  <*>u  in  Eq.  (14)  are  zero  since  u  extends  only  over  (-1/2,  1/2).  Then  K 
equals  2P,  that  is,  regardless  of  the  window,  the  EDF  increases  linearly 
with  the  number  of  pieces  P  until  overlap  occurs  (see  Eq.  (13)).  As  P  in¬ 
creases  somewhat  beyond  T/L,  K  continues  to  increase,  although  at  a  slower 
rate,  because  the  overlapped  pieces  are  progressively  more  statistically  de¬ 
pendent.  Of  importance  in  the  behavior  of  K  are  the  rate  of  increase  of  K 
with  P,  and  the  maximum  value  of  K  attainable  through  the  choice  of  P. 

LIMITING  VALUE  OF  EQUIVALENT  NUMBER  OF  DEGREES  OF  FREEDOM 

As  the  number  of  pieces  P  tends  to  infinity,  the  overlap  approaches  100 
percent,  and  the  denominator  of  Eq.  (14;  approaches  an  integral,  yielding 


•> 
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This  limiting  value  depends  only  on  T/L  and  the  correlation  <*>u  of  the  shape 
u  of  the  window.  It  is  finite  because  the  overlapped  pieces  are  statistically 
dependent. 

For  large  values  of  T/L  (ratio  of  record  length  to  segment  length),  an 
alternate  form  of  Eq.  (15)  is  very  illuminating.  If  Eq.  (10)  is  utilized,  the 
denominator  approaches 


f,  Vl1  A  VLl>  _i_f 

J  «.,(»)  'J  v0)  ~LJ 


<t>  \l) 

w 

v°> 


(16) 


but  from  Eqs.  (7)  and  (4), 


and 


Riving 


<*>w(0)  =/dt  |w(t)|2  =/df  |W(f)j 2 , 
<*>w(t)  -  j~di  exp  (i  2  n  ft)  |\V(f)|2  , 

f  dt  |Vt)|2=/df|w<f)|4 


by  Parseval's  Theorem.  Combining  Eqs.  (15)  through  (19),  we  obtain 

.12 


K  S  2  (T  -  LH 

<I> 


for  — » 1  . 

i-i 


(17) 

(18) 

(19) 


[/-f  tw(f)|2J 

y*dfiw(f)! 

If  wc  define  the  statistical  bandwidth  [5.  p.  263]  of  spectral  window  |  W|^  as 

[/dilwwl2]  _cst 
st  y*df!\v(oi4  L 


B  . 


then 


K  (T  -  L)  B  =  2  (~  -  l)  C 
<*>  st  '  L  /  st 


for  ~  >:>  1  • 

L 


(22) 
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The  constant  Cst  is  dimensionless  and  of  the  order  of  unity;  it  deponds  only  on 
the  shape  of  the  window: 


C 


st 


where  U  is  the  Fourier  transform  of  u  (see  also  Eq.  (9)). 


(23) 


The  first  form  in  Eq.  (22)  for  Ka>  indicates  that  if  windows  are  compared 
on  the  basis  of  equal  statistical  bandwidths  (by  appropriate  choice  of  segment 
length  L  for  each  data  window),  then  all  windows  have  the  same  value  of  KU1 
for  large  T/L;  that  is,  all  windows  have  the  same  variance-reduction  capabil¬ 
ity,  when  compared  on  the  basis  of  equal  statistical  bandwidths,  if  the  avail¬ 
able  record  length  is  much  larger  than  the  segment  length.  For  other  measures 
of  bandwidth,  such  as  the  half -power  bandwidth  B,  we  are  led  to  anticipate 
this  same  result.  In  a  later  section,  we  will  demonstrate  this  quantitatively  not 
only  for  T/L  »  1  but  for  small  values  of  this  ratio  as  well,  and  for  finite 
values  of  P,  the  number  of  pieces. 


It  is  of  interest  to  compare  Eq.  (20)  with  the  results  of  Blackman  and  Tukey 
[2,  Secs.  B6-B8]  for  spectral  estimation  via  the  "indirect"  correlation  func¬ 
tion  approach.  Their  EDF  at  frequency  fj  is  given  approximately  by 


2T 


r  •<ja 

L/df  [Qi(f+f1)  +Qi(f-f1)]  G(f)J2 
df  [Q.ff-K^  +Qi(f-f1)]2  G2(f) 


(24) 


for  long  records,  where  is  their  spectral  window.  For  frequencies 
greater  than  the  width  uf  spectral  window  Qj,  and  assuming  that  Qi  »s  narrow 
compared  with  the  finest  detail  in  true  spectrum  G,  Eq.  (24)  becomes  approxi¬ 
mately 
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In  order  to  relate  the  EDF  in  Eq.  (25)  to  the  one  used  here,  we  note  that  the 
._an  value  of  Blackman  and  Tukcv's  spectral  estimate  is  gi  ven  by  the  convolu- 
on  of  with  G.  We  then  identify  Qj  with  |wp,  obtaining  for  the  EDF  in 
Cq.  (25). 

[/W]* 

2T  — ^ - ,  (2G) 

I  df|W(f)|4 

-a> 

which  is  in  agreement  with  Eq.  (20)  for  long  records.  Thus,  under  the  assump¬ 
tions  given  above,  the  same  limiting  value  of  EDF  is  realized  by  both  the 
"indirect"  correlation  approach  and  the  present  "direct"  FFT  approach,  that  is, 
both  methods  are  capable  of  the  same  statistical  stability  if  the  proper  overlap 
is  used  in  the  FFT  approach. 


DATA  WINDOWS  AND  CHARACTERISTICS 

Four  data  windows  will  be  considered  here.  They  are  all  continuous;  how¬ 
ever,  they  have  differing  degrees  of  continuity  in  their  derivatives,  leading  to 
different  rates  of  decay  of  their  spectral  windows  for  large  frequencies.  The 
triangular  data  window,  made  up  of  two  straight-line  segments,  has  a  discon¬ 
tinuous  first  derivative.  The  cosine  data  window  (Hanr.ing)  has  a  discontinuous 
second  derivative.  The  quadratic  data  window  is  made  up  of  segments  of 
quadratic  curves  so  chosen  that  the  first  derivative  is  everywhere  continuous, 
but  the  second  derivative  is  discontinuous;  thus,  the  quadratic  data  window  has 
behavior  similar  to  that  of  the  cosine  data  window.  The  cubic  data  window  is 
made  up  of  segments  of  cubic  curves  so  chosen  that  the  second  derivative  is 
everywhere  continuous,  but  the  third  derivative  is  discontinuous.  This  se¬ 
quence  of  windows  will  have  progressively  better  high-frequency  decay  and,  as 
will  be  seen  shortly,  better  side  lubes  at  low  frequencies.  (Hamming  weighting 
is  not  considered  because  its  discontinuous  data  window  yields  a  slowly  decay¬ 
ing  spectral  window  for  large  frequencies . )  Computation  of  the  quadratic 
and  cubic  data  windows  is  easier  and  faster  than  computation  of  the  cosine  data 
window  in  the  time  domain.*  Their  computational  advantage  and  better  side- 
lobe  behavior,  make  them  attractive  candidates  for  spectral  analysis. 

*ln  the  frequency  domain,  the  cosine  data  window  is  equivalent  to  convolu¬ 
tion  with  the  sequence  -1/4,  1/2,  -1/4,  which  is  easily  implemented. 
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The  four  windows  arc  detailed  below  in  Eqs.  (27)  through  (30).  They  have 
been  normalized  in  such  a  way  that  U(0)  =  1.  From  Eq.  (9)  end  Fig.  1,  notice 
that  u(t)  =  0  for  it  I  >  1/2.  (The  expressions  for  the  correlations  <f>u  of  the 
windows  are  collected  in  Appendix  B,  these  correlations  are  necessary  for  the 
evaluation  of  Eq.  (14).)  In  the  following,  sinc(x) s  sin(rx)/(TX). 


Triangular  Data  Window 


u(t>  =  2  (1  -  2ltl),  ltl<~ 


U  (f)  =  sine  (f/2)  , 


(2  7) 


Cosine  Data  Window 


u(t)  =  1  +  COS(2rrt).  1 1 1  <  — 


U(f)  =  ^  . 

i  -  r 


(28) 


Quadratic  Data  Window 


’J(t)  = 


,t.Sl 

T(‘-2't')2’  2  ) 


U(f)  =  sine  (f/3)  . 


(29) 


Cubic  Data  Window 


u(t)  = 


~(l  -  24t2  +  48lt|3),  It l  <  ~ 


U(f)  -  sine  (f/4)  . 


(30) 


The  function  \V(f)  is  available  from  the  normalized  function  U(f)  accord¬ 
ing  to 

W  (f )  =  L  U  (Lf)  ,  (31) 

upon  Fourier  transformation  ol  Lq.  (9).  We  define  the  half-power  bandwidth 
of  spectral  window  |W|2  as  the  frequency  range  over  which  |W|2  is  greater 
than  half  of  its  peak  value. 


|w(*|b)|2=!|w<o)|2.  t^> 

In  addition  to  the  statistical -bandwidth  constant  Cst  delined  in  Eq.  (21),  we 
define  a  half-power-bandwidth  constant  C  according  to 

B  =  £  •  (33) 

Numerical  values  for  both  of  these  dimensionless  constants  for  the  four  win¬ 
dows  are  given  in  Table  1. 

The  bandwidth  constants  are  larger  for  the  "smoother"  data  windows;  thus, 
their  bandwidths  arc  larger  for  a  given  segment  length  L  (Fig.  1).  Alternatively, 
If  the  bandwidths  arc  to  be  kept  equal  for  the  four  windows,  the  segment  lengths 
must  be  larger  for  the  smoother  data  windows. 

In  the  bottom  row  of  Table  1,  the  ratio  of  the  statistical-bandwidth  constant 
to  the  half-power-bandwidth  constant  is  found  to  be  relatively  constant  for  the 
four  windows  considered.  Thus,  the  statistical  stability  of  the  spectral  esti- 


12 


Table  1 


BANDWIDTH  CONSTANTS 


Data  Window 

Triangular 

Cosine 

Quadratic 

Cubic 

cst 

1.854 

2.070 

2.304 

2.  686 

C 

1.27G 

1.441 

1,572 

1.820 

Cst'C 

1.454 

1.443 

1.4GG 

1.476 

mates  can  be  discussed  in  terms  of  either  bandwidth  without  tear  of  changing 
significantly  the  quantitative  aspects.  For  example,  Eq.  (22)  becomes 

2.9(T-L)B,  for  £  »  1,  (34) 

in  terms  of  half-power  bandwidth  B.  where  Eqs.  (21)  and  (33)  and  Table  1 
have  been  employed.  A  more  precise  relation  than  Eq.  (34)  will  be  given  for 
K  in  the  next  section,  where  the  number  of  pieces  P  will  be  finite. 

Half  of  each  symmetric  spectral  window'  is  plotted  in  dB  versus  f/B  in 
Figs.  3  through  G.  Here 

dBElO  log^lul2,  (35) 

since  the  power-density  spectrum  G  is  seen  through  a  window  proportional 
to  |U|2.  All  the  plots  go  through  -3. 01  dB  at  f/B  =  1/2,  since  B  is  the  half¬ 
power  bandwidth.  The  slow  spectral  decay  of  the  triangular  data  window  and 
the  fast  spectral  decay  of  the  cubic  data  window  are  evident.  The  cosine  and 
quadratic  data  windows  exhibit  intermediate  behavior.  The  first  three  side 
lobes  of  the  spectral  windows  are  given  in  Table  2.  where  it  is  seen  that  the 
quadratic  window  offers  an  8.3-dB  improvement  relative  to  the  cosine  window' 
in  the  size  of  the  first  side  lobe,  and  the  cubic  window  yields  an  additional 
13. 3-dB  improvement.  . 


Spectral  Window  for  Cosine  Data  Window 


Fig.  5.  Spectral  Window  for  Quadratic  Data  Window 


Fig.  6.  Spectral  Window  for  Cubic  Data  Window 


Table  2 


FIRST  THREE  SIDE  LOBES  OF 
SPECTRAL  WINDOWS 


Data  Window 

Side  Lobes  (dB) 

Triangular 

-26.5 

-35.7 

-41.6 

Cosine 

-31.5 

-41.5 

-48.5 

Quadratic 

-39.  8 

-53.  5 

-62.4 

Cubic 

-53.  1 

-71.3 

-83.2 

RESULTS 


The  general  expression  for  EDF  is  given  in  Eq.  (14).  Wc  eliminate  the 
segment  length  L  in  this  expression  in  favor  of  the  half-power  bandwidth  B 
by  using  Eq.  (33)  to  obtain  the  dependence  on  the  fundamental  parameter  BT 
(see  Introduction).  It  follows  that 


P-1 

1  - 

*  (khJ/c-‘) 

2 

E  ( 

u'  P-1  / 

1  p/ 

4  (6) 

.  s 

1 

It 

u 

(36) 


For  a  particular  window  u  and  value  of  BT,  K  is  computed  versus  F.  A 
sample  tabulation  for  the  cosine  data  window  and  BT  -  8  is  given  in  Table  3. 
The  column  headed  "Fractional  Overlap"  is  a  measure  of  how  much  the  individual 
data  windows  overlap  in  the  spectral  processing  technique  depicted  in  Fig.  1  and 
is  given  by 
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Table  3 


EQUIVALENT  DEGREES  OF  FREEDOM 
VERSUS  NUMBER  OF  PIECES, 

BT  =  8,  COSINE  DATA  WINDOW 


p 

K 

Fractional 

Overlap 

2 

4.00 

.00 

3 

6.00 

.00 

4 

8.00 

.00 

5 

10.00 

.00 

G 

12.00 

.09 

7 

14.00 

.24 

8 

15.96 

.  35 

9 

17.74 

.43 

10 

19. 13 

.49 

11 

20.03 

.54 

12 

20.50 

.59 

13 

20.  09 

.62 

14 

20.  73 

.65 

15 

20.71 

.67 

16 

20.  66 

.70 

17 

20.  61 

.72 

18 

20.56 

.73 

19 

20.52 

.75 

20 

20.48 

.76 

30 

20. 19 

.84 

40 

20.  04 

.  88 

50 

19.  95 

.91 

100 

19.77 

.  95 

200 

IS.  67 

.98 

(37) 


L-S  P-T/L  P-BT/C 
L  P-1  P-1 


when  non-negative.  Equations  (13)  and  (33)  were  employed  in  Eq.  (37). 


Several  points  In  Table  3  are  worth  noting.  For  no  overlap,  the  EDF 
increases  linearly  with  P,  and  even  as  overlap  begins  to  occur,  the  rate 
of  increase  of  EDF  remains  the  same.  Thus,  249  overlap  still  yields  the 
maximum  possible  EDF  for  P  =  7.  However,  as  P  and  the  overlap  increase 
further,  the  EDF  increases  more  slowly,  eventually  reaching  a  maximum,* 
after  which  it  decreases  slightly  for  further  increases  in  P.  A  point  of 
diminishing  returns  is  reached  somewhere  near  P  -  12,  where  989c  of  the 
maximum  (max)  EDF  is  realized.  The  extra  computational  effort  in  spectral 
analysis  for  P  >  12  is  not  worth  the  return  in  stability.  The  fractional  over¬ 
lap  for  989  of  max  EDF  is  .59;  the  EDF  is  then  20.5,  whereas  it  was  only 
10.0  for  the  last  non-overlapped  example.  The  case  for  overlapped  processing 
is  well  demonstrated  by  Table  3. 


Similar  ..-_ulto  for  all  four  windows  and  BT  -  2,  4,  8,  16,  32,  and  64 
are  condensed  in  Table  4,  which  gives  the  required  number  of  pieces  and 
the  corresponding  fractional  overlap  for  a  specified  fraction  of  the  max  EDF. 
For  example,  in  Table  4B  for  the  cosine  data  window,  in  order  to  realize 
a  specified  fraction  of  989  of  the  max  EDF  at  BT  =  8,  the  number  of  pieces 
required  is  12,  and  the  corresponding  fractional  overlap  is  .59.  The  bottom 
row  of  each  data-window  table  gives  the  max  EDF  for  the  corresponding 
value  of  BT.  Also  provided  is  an  equation  for  max  EDF,  which  was  empiri¬ 
cally  determined  to  fit  through  the  numerical  values  obtained. 


Several  striking  invariances  are  apparent  upon  inspection  of  Table  4.  First, 
for  a  given  record  length  T  and  desired  frequency  resolution  B,  the  max 


♦The  existence  of  a  finite  value  of  P  for  maximum  EDF  is  similar  to  the 
situation  cited  in  Reference  11. 
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EDF  is  virtually  independent*  of  the  window  employed;  that  is,  all  the  windows 
considered  have  the  same  variance-reduction  capability  in  spectral  estimation 
when  compared  under  the  same  frequency- resolution  constraint.  A  simple, 
approximate  rule  of  thumb  for  all  four  windows  is  given  by 

max  EDF  w3  (BT  -  1).  (38) 

Recall  that  B  is  the  half-power  bandwidth  of  the  spectral  window. 

For  a  given  window  and  specified  fraction  of  the  max  EDF,  it  will  be  ob¬ 
served  from  Table  4  that  for  BT  >  4,  the  required  fractional  overlap  is 
approximately  constant  (independent  of  BT).  Therefore,  in  the  last  column  of 
Table  4  is  entered  a  representative  or  "average”  fractional  overlap  required 
for  the  specified  fraction  of  max  EDF  entered  in  the  first  column.  (This  simple 
rule  does  not  hold  well  when  100Cc  of  the  max  EDF  is  required;  accordingly,  no 
value  is  entered  tor  this  case.  ) 


For  a  given  value  of  BT,  the  number  of  pieces  required  to  realize  .9 9  (or 
less)  of  max  FDF  is  v  irtually  independent  of  the  particular  window  employed  in 
spectral  analysis.  This  independence  is  very  important;  it  says  that  all  four 
windows  require  the  same  number  of  FFTs  in  order  to  realize  the  same  EDF 
and  that  selection  among  the  •  therefore,  can  not  be  based  upon  the 

number  of  FFTs  required,  1  based  upon  some  other  consideration 

such  as  side-lobe  level  or  s  (to  be  discussed).  An  approximate  rule 

of  thumb  for  the  number  oi  p.  ,«  ired  is  given  by 

Number  of  pieces  required  ■'■1.75  BT  for  .  99  of  max  EDF.  (39) 

For  large  BT  products,  the  number  of  pieces  required  to  realize  max  EDF 
is  significantly  larger  than  the  number  required  to  realize  99 .  of  max  EDF; 
thus,  this  large  amount  of  additional  processing  yields  insignificant  improve¬ 
ment  and  is  to  be  avoided. 

•Very  small  values  of  BT  arc  an  exception;  these  are  of  little  practical 
interest  however. 


Tabic  1C.  Quadratic  Data  Window 


Table  5 


REQUIRED  FRACTIONAL  OVERLAP  FOR  .99  max  EDF 


Data  Window 

Triangular 

— 

Cosine 

— 

Quadratic 

Cubic 

Fractional 

Overlap 

.  56 

.61 

.65 

.70 

The  required  fractional  overlap  is  greater  for  the  hetter-side-Iobe  win¬ 
dows.  Thus,  for  example,  to  realize  .99  of  max  EDF,  we  have  (virtually  inde¬ 
pendent  of  the  BT  product)  the  values  listed  in  Table  5.  Notice  that  rather 
large  overlaps  are  required  for  some  windows. 

The  difficulty  of  realizing  general  fractional  overlaps,  such  as  .56,  raises 
the  question  as  to  what  fraction  of  max  EDF  L  attainable  if  one  restricts  over¬ 
laps  to  a  few  easily  realized  overlaps  such  as  .50  and  .625.  This  question  Is 
answered  in  Tables  6  and  7.  *  Table  6  indicates  that  the  cosine  data  window  at 
50  overlap  (a  popular  case)  yields  92'/  of  the  max  EDF.  However,  the  cubic 
data  window  realizes  only  7 59c  of  its  potential  at  5(Tr  overlap.  Table  7  shows 
that  when  the  fractional  overlap  is  increased  to  5/8,  the  cosine  and  quadratic 
data  windows  realize  virtually  their  ultimate  capability.  If  the  overlap  is  in¬ 
creased  to  7 5'  (  ,  the  cubic  data  window  then  realizes  its  max  EDF. 


Table  6 

ATTAINABLE  FRACTION  OF  max  EDF 
AT  .  50  FRACTIONAL  OVERLAP 


Data  Window 

Triangular 

Cosine 

Quadratic 

Cubic 

Fraction  of 
max  EDF 

.96 

.92 

.85 

.75 

‘These  values  are  not  attainable  from  Table  4  but  come  from  the  complete 
tabular  results,  of  which  Table  3  is  one  example. 
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Table  7 


ATTAINABLE  FRACTION  OF  MAX  EDF 
AT  .625  FRACTIONAL  OVERLAP 


Data  Window 

Cosine 

Quadratic 

Cubic 

Fraction  of 
max  EDF 

1.00 

.  98 

.93 

Thus  far,  no  trade-off  has  been  necessary  to  realize  the  better  side  lobes 
of  the  smoother  data  windows;  that  Is,  by  proper  choice  of  overlap,  equal  statis¬ 
tical  stability  is  attainable,  and  an  equal  number  of  FETs  is  required,  for  all 
four  windows.  However,  there  is  one  trade-off  that  enters  as  follows;  if  the 
original  record  length  T  is  composed  of  samples  of  a  process  at  increments* 
At,  more  samples  are  contained  in  the  segment  length  L  for  the  better-side- 
lobe  windows;  that  is,  th™  number  of  samples  in  interval  L  is 


N 

s 


_L _ C__ 

At  ~  B  At 


(4  0) 


employing  Eq.  (33).  Ns  is  directly  proportional  to  half-power-bandwidth 
constant  C  for  a  specified  sampling  increment  At  and  resolution  B.  Thus, 
using  Table  1,  the  cosine  data  window  requires  1. 13  times  as  many  samples  as 
the  triangular  data  window  requires.  The  corresponding  ratios  for  the  quad¬ 
ratic  and  cubic  data  windows  are  1.23  and  1.43,  respectively.  Thus,  better 
side  lobes  in  spectral  analysis  can  be  realized  at  the  expense  of  larger-size 
FFTs,  rather  than  at  the  expense  of  statistical  stability  or  number  of  FFTs. 
These  comments  hold  for  equal  half-power  bandwidths  of  the  windows. 

For  a  specified  sampling  increment  At,  desired  resolution  B,  and 
particular  w  indow  ,  Eq.  (40)  will  generally  not  be  a  power  of  2.  Since  FFTs  run 


‘The  sampling  increment  At  must  be  chosen  small  enough  to  avoid  alias¬ 
ing;  this  is  the  only  area  where  the  bandwidth  of  the  process  comes  into  con¬ 
sideration. 
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faster  when  conducted  at  powers  of  2,  it  is  recommended  that  the  desired 
resolution  B  be  changed  somewhat  (increased  or  decreased)  so  as  to  make  Ns 
a  power  of  2.  This  is  generally  a  tolerable  situation  since  B  is  often  a 
"guesstimate"  in  the  first  place. 


OPTIMUM  WEIGHTING  OF  INDIVIDUAL  SPECTRAL  ESTIMATES 

In  Eq.  (3),  the  p-th  estimate  |Yp(t)j2  of  the  power-density  spectrum  was 
weighted  equally  with  all  other  estimates.  In  this  section,  .  e  consider  whether 
unequal  weighting  will  yield  additional  worthwhile  variance  reduction.  Inas¬ 
much  as  the  edge  pieces  in  Fig.  1  are  weighted  only  once  by  a  data  window, 
whereas  the  interior  pieces  are  weighted  mote  than  once,  perhaps  heavier 
weighting  of  the  edge  pieces  will  yield  additional  stability.  The  power-density 
estimate  is  formed  as 


P 

to  *2  “plv'l2-  <41> 

p=l  F  '  F  ' 

The  derivation  of  the  EDF  of  this  estimate  is  given  in  Appendix  C,  which  also 
presents  the  optimization  of  the  EDF  by  choice  of  weights  for  a  given  P, 
record  length,  and  window.  A  summary  of  the  numerical  results  is  given  in 
Table  8,  where  P  is  varied  up  to  64.  The  largest  value  of  EDF  attained  over 
that  range  of  P  is  quoted  In  Table  8,  except  for  BT  =  16  where,  with  the  ex¬ 
ception  of  the  triangular  data  window,  P  =  64  was  not  yet  great  enough  tc 
reach  the  max  EDF  by  weighting. 


Table  8 

OPTIMUM  EDF  VALUER 


Data  Window 

BT 

2 

4 

8 

10 

Triangular 

4.37 

10.  12 

21.65 

44.72 

Cosine 

4.50 

10. 15 

21.44 

>43.87 

Quadratic 

4.07 

9.80 

21.29 

>44. 18 

Cubic 

3.82 

9.  52 

20.  94 

> 43.72 

2G 


A  comparison  ot  Table  8  with  the  max  EDF  values  of  Fig.  4  reveals  that 
very  little  is  to  be  gained  bv  optimum  weighting,  except  for  small  values  of  BT. 
However,  small  values  of  BT  are  not  of  great  practical  interest  because  the 
estimates  are  very  unstable  statistically.  Also,  the  number  of  pieces  P 
required  to  realize  the  optimum  EDF  is  rather  large;  for  example,  in  order 
to  gain  an  improvement  in  EDF  of  0.5  over  the  max  EDF  in  Table  4,  19  pieces 
are  required  for  the  cosine  data  window,  20  pieces  for  the  quadratic  data  win¬ 
dow,  and  15  pieces  for  the  cubic  data  window.  Moreover,  the  optimum  weights 
are  found  to  alternate  in  sign  for  some  cases,  causing  a  loss  in  significance. 


DISCUSSION 

This  investigation  of  four  good  data  windows  indicates  that,  there  is  no  best 
window  for  spectral  estimation.  Rather,  there  is  a  trade-off  to  be  made  when 
choosing  a  w  indow  :  the  better-side-lobe  w  indow  s  require  larger-size  FFTs. 
When  the  proper  overlap  is  used  for  each  data  window,  the  selection  of  u  in- 
dows  can  not  be  made  on  the  basis  of  statistical  stability  or  the  number  of  FFTs 
required. 

The  quadratic  and  cubic  data  windows  are  simpler  and  quicker  to  compute 
than  the  cosine  data  window  in  the  time  domain  (but  not  in  the  frequency  do¬ 
main  as  regards  their  effects).  In  addition,  the  quadratic  and  cubic  windows 
have  better  side-lobe  behavior  and,  therefore,  merit  serious  consideration  for 
spectral  analysis.  However,  they  require  larger-size  (but  not  move)  FFTs 
than  does  the  cosine  data  w  indow  . 

The  reason  that  the  better-side-lobe  windows  do  not  require  more  FFTs 
than  do  the  other  windows  is  as  follows.  For  a  fixed  half-power  bandwidth  B, 
the  better-side-lobe  windows  require  larger  segment  lengths  L;  however,  the 
corresponding  data  windows  tend  to  be  more  peaked  near  the  center  of  the  seg¬ 
ment  length.  In  order  to  utilize  a  given  record  length  for  maximum  statistical 
stability,  these  data  windows  must,  therefore,  overlap  for  a  greater  percentage 
of  the  segment  length.  It  turns  out  that  the  increased  segment  length  and 
increased  overlap  almost  exactly  compensate  each  other,  so  that  a  constant 
number  of  FFTs  is  required  regardless  of  the  window  selection. 

This  report  has  concentrated  on  the  variance  of  the  spectral  estimates.  In 
Appendix  D,  it  is  shown  that  the  covariance  of  spectral  estimates  at  two  differ¬ 
ent  frequencies  is  always  positive  but  is  essentially  zero  when  the  frequencies 
differ  by  more  than  the  w  idth  of  the  spectral  w  indow  .  Thus,  spectral  estimates 
at  frequencies  farther  apart  than  B  are  statistically  linearly  independent  of 
each  other. 


27  '28 

REVERSE  BLANK 


Appendix  A 


DERIVATION  OF  MEAN  AND  VARIANCE 


From  Eqs.  (2)  and  (3),  the  spectral  estimate  at  frequency  f  is 
G(f)  =  Jjf  iu  dv  exp(-i2*f(u-v))  x(u)x*(v)  • 

w[u  -i-  -  (p-l)S]  w*[v  -  y  -  (p-l)S]  . 


(A  —  I) 


(For  completeness,  process  x  and  window  w  are  allowed  to  be  complex.) 
The  mean  value  of  u(f)  is  obtained  by  ensemble-averaging  Eq.  (A-l)  over  the 
possible  realizations  of  process  x.  Expressing  the  correlation  of  x  as  a 
Fourier  transform  of  spectrum  G,  we  find  the  average  value 

P 

EjCl(f)}  =-p-^  jjdudv  exp(-i2»f(u-v)J J'do  exp(i2ir(u-v>)  G(y)» 
w[u  -  y  -  (p-l)S]  w*  [v  -  ~  -  (p-l)s] 


-/*  G(*)l\V(f-,)|2  \f'h  G(f- v)  |\V(,)|2  , 


(A-2) 


where  we  have  utilized  Eq.  (4).  If  spectral  window  ItVl2  is  narrower  than  the 
finest  detail  in  the  true  spectrum  G,  Eq.  (A-z)  becomes 


E|G(f)|  =  G(t) y* dt'  |W(^)|2  . 


(A -3) 


Equation  (A-2)  is  not  limited  to  Gaussian  processes  but  is,  in  fact,  true  for  any 
stationary  process. 

In  order  to  evaluate  the  variance  of  the  spectral  estimate,  we  start  with 


Prororiino  nacre  hlanir 

•  ivuwiilii^  uiy lilt 
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\V[u.i.(p_1)S]  w*[v-j-(p-l)S]  w[r-j-(q-l)S]  w*[s -  (q-l)S]  . 

(A -4) 

In  order  to  simplify  this  expression,  we  must  be  able  to  evaluate  the  fourth- 
order  average.  If  x  is  a  real  Gaussian  random  process,  the  average  in  Eq. 
(A-4)  becomes 


Il(u-v)  R(r-s)  +  R(u-r)  R(v-s)  +  R(u-s)  R(r-v)  ,  (A-3) 

where  R  is  the  correlation  of  process  x.  (If  x  is  a  complex  envelope  of  a 
Gaussian  process,  the  middle  term  in  Eq.  (A-5)  is  absent  [12]. )  When  we 
express  correlation  R  as  a  Fourier  transform  of  spectrum  G  and  substitute 
Eq.  (A-5)  in  Eq.  (A-4),  we  obtain 


+  exp(i2jTAi(u-s)+i27ri-(r--v)]] 


tit* 


di'  G(m)  G(»)  [l\V(f-M)|  I W (f  —  e)|^  +  cxp(i27r(M"  '•■)(p-q)S)* 


W(f-M)  W*(f+  v)  W’(ft-M)  W*(f-K)  +  cxp(i2ir(M- •')(p-q)S)  |\V(f-M)|“  |\V(f-t')|^]  , 

(A -6) 

using  Eq.  (4).  The  first  term  in  Eq.  (A-G)  is  recognized  from  Eq.  (A-2)  as  the 
square  of  the  mean  of  u.  Therefore, 


:i0 


Var{6(f)|  = 


■ir-m 


d M  G(u)  exp(i2xM(p-q)S)  \V(f->d  W(f+A») 


P=i  q; 


■htt  fdn  G(/i)  exp(i2fAi(p-q)S)  t\V(f-M)|2|  .  (A-?) 

^1  a?l  J  ' 


P  p=l  q=l 


Now  if  f,  the  frequency  of  interest,  is  greater  than  the  bandwidth  of  the  spec¬ 
tral  window  |\v|~  (i.e. ,  a  couple  of  resolution  cells  away  from  the  origin), 
then  W(1-m)  and  \V(f+a)  do  not  overlap  significantly.  Letting  B  be  the  half¬ 
power  bandwidth  of  the  spectral  window  1W|2,  we  therefore  have  the  excellent 
approximation*  for  f  >B, 


VarjG(f)}  = 


^  ^  ^JdnG(n)  exp(i2ir^(p-q)S)  l\V(f-M)l2 


(A-8) 


(For  1  =  0,  the  two  terms  in  Eq.  (A-7)  are  equal  if  data  window  w  is  real,  in 
which  case  Var[O(0)}  is  double  that  given  by  Eq.  (A-8)  at  f=0.)  Making  he 
change  of  variable  k  -  p-q  in  Eq.  (A-8),  wc  obtain 


P-1  \  /* 

VarjG'(f)}  S-i  V  (l  -  ~)  /d*.  exp(i2^kb)  G(m)  |W(f-M)|2 

^-(P-ir  V,\J 


.  (A -9) 


But  if  the  bandwidth.  B  of  spectra]  window  |\V|^  is  narrower  than  the  lines, 
detail  in  spectrum  G,  the  integral  on  n  in  £q.  (A -9)  can  be  approximated  by 

C(f)  /da  exp(i2*-MkS)  I W (f — m)I  2 


=  G(f)  exp(i2iifkS)  <t>  (kS)  , 
w 


(A-10) 


where  we  have  utilized  Lqs.  (4)  and  (7).  Then  Eq.  ( 


'  K  f\\  l _ _ _ 


*\Vhen  x  is  a  complex  envelope,  Eq.  (A-8)  is  exact;  see  comment  under 
Eq.  (A -5). 
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(A -11) 


Vari%)}2G2(f)|  £  f1  *  7^w • 

k  =  -(P-1) 

Ihis  equation  is  similar  to  that  given  in  [9,  p.  ?l]  . 

Equations  (A -3)  and  (A-ll)  are  the  main  products  of  this  appendix. 
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Appendix  B 


CORRELATIONS  OF  DATA  WINDOWS 


The  correlation  of  the  normalized  data  window  u  is  given  by  Eq.  (11)  as 


Vr) 


■/* 


=  /  dt  u(t)  u*(t-T> 


(B-l) 


This  quantity  is  required  for  the  calculation  of  EDF,  K,  in  Eq.  (14).  Since 
u(t)  =  0  for  iti  >1/2,  then  4>u(r)  =  0  for  m  >  1.  Thus,  for  the: 


Triangular  Data  Window 


2  3  1 

<t>  (r)  |  1  -  6r  +  6 iti  ,  |r|  <  - 

<t>  (0)  i  ■>  i 

U  I  2  (1  -  |T|  )  ,  -  <  |T|<  1 


(B-2) 


Cosine  Data  Window 


^u1^  2  f  1  11 

=  “  (1  -  m  )  jl  +  JCOS(2irr)  +—  Sill(2irir|),  i t i  <  1 


<j»u(0)  3 


(B-3) 


Quadratic  Data  Window 


<t> 


<t>  (0) 
u 


81  5  2 

Q  (r)  -  — — -  It i)  °  =  Q  (r),  ~  <  l7l  <  — 

11  \3  /  2y  3  3 


(B-4) 
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Cubic  Data  Window 


Vr> 
*  (0) 


„  ,  28672/1  \  7  1  ^  ,  1 

C  <  r)  + - (—  -  |T|)  =;  C  (r),  -  <rn<  — 

2'  '  151  \ 2  7  3l  ''  4  _  2 


57344/1 


l03(,,"isi 


ITI  <  - 


(B-5) 


The  forms  for  the  correlation  of  the  quadratic  and  cubic  data  window  a  are 
compact  and  very  useful  for  computer  programming. 
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Appendix  C 


where 


0  <kS) 
w 


V0) 


(C-6) 


Partially  differentiating  K  with  respect  to  wj1  [l3,  Appendix]  and  noting  that 
the  absolute  scale  of  weights  jwpj  does  not  enter  in  R,  wc  sec  that  the 
optimum  weights  must  satisfy  the  equation 


E 


M.  w  =  1, 
k  j  +  k 


1  <  j  <  P  . 


(C  -  7) 


If  we  define  the  matrices 


,\4 


[m 


1 

m-nj 


1  <  m ,  n  <  P  , 


1T  =  [)  "...  l],  wT=  [Wl  w2  ...  wp]  ,  (C-U| 


then  the  optimum  weights  are 


w  -  M~  ^  1  > 


and  the  optimum  EDF  is 


K  =  2  lTM'll  . 


(C-») 


(C-10) 


:u; 


gr~^r--ryj 


Appendix  D 


COVARIANCE  OF  SPECTRAL  ESTIMATES 


A  generalization  of  the  technique  in  Appendix  A  leads  to  the  following 
expression  fo*  the  covariance  between  spectral  estimates  at  frequencies  fj  and 

f2  : 


Cov 


1%).  6(f2>l  =  Z)  Z|  |  fdti  G(M)  W(fr 

p  p 

z 

q= 


M)  W(f2+^i)  exp(i2?r/i(p-q)S) 


,  P  p  i  r 


dM  G(m)  W(f  -m)  W*(f  ->j)exp(i2,.H(p- 


q)h|  . 

(D-l) 


Now  if  fj+f2  is  greater  than  the  bandwidth  of  the  spectral  window,  W(f^-ji) 
and  W(f2+M)  are  essentially  non-overlapping.  Then, 


Cov  =  ~ 


P-l 

2 

k=-(P-i) 


<-W 


dn  G(m)  Wffj-n)  W*<f  -/i)  exp(i2rukS) 


-(D-2) 


This  quantity  is  always  positive.  However,  it  is  very  small  when  | f2~f x ) > B 
because  W(fj-ji)  and  W(f2->d  do  not  overlap  then. 

When  true  spectrum  G  varies  but  slightly  over  a  frequency  range  B, 


P-1 

Cov  SGtfj)  G(f2»  £  (l  -■—) 

k=-(P-l) 


|x(kS ,  f2- f^ 


(0-3) 


if  7 


where 


X(r 


f’5/' 

-/ 


d^i  exp(i2»r>ir)  W(x+f)W*(M) 
dt  exp  (-i2»rft)  w(t)  w*(t-r) 


(D-4) 


is  the  ambiguity  function  of  window  w  .  Again,  if  |  fi-f2  |  >  B  ,  Eq.  p-3)  is 
essentially  zero,  as  shown  by  the  first  form  in  Eq.  (D-4). 


As  an  example,  for  the  cosine  window  and  01.  overlap,  the  covariance 
coefficient  (ratio  of  Cov  to  the  square-root  of  the  product  of  variances)  is 

9  ’  36  *  °.  °.  •  •  •  •  for  |f2~fi|  =  0,L’L’L*L . respectively. 

2 

Thus  spectral  estimates  —  Hz  apart  are  essentially  uncorrelated.  For  501 

overlap  (and  large  P),  the  corresponding  covariance  coefficients  are  slightly 
larger,  being  1,  ,195,  .068,  .005,  0 . 
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